Phase diagrams of a Spin-1 Ising Superlattice 

Hamid EZ-ZAHRAOUY^^and Ahmed KASSOU-OU-ALI^ 

April 29, 2004 

(a) Faculte des Sciences, Departement de Physique, Laboratoire de Magnetisme et Physique 
des Hautes Energies, B.P. 1014, Rabat, Morocco, 
(b) Faculte des Sciences, Departement de Physique, Laboratoire de Physique Theorique, 

B.P. 1014, Rabat, Morocco. 

Abstract 

The three-dimensional spin-1 Ising superlattice consisting of two different 
ferromagnetic materials with two different crystal fields Ai and A2 is considered 
in the mean field approximation. The phase diagrams are considered in the 
(t,c?2) plane for different ranges of variation of d\ (t=T/J, cfi/J and cfo/J are the 
reduced temperature and crystal fields respectively). The phase diagrams exhibit 
a variety of multicritical points and reentrant and double reentrant behaviours. 
They are found to depend qualitatively and/or quantitatively on the thicknesses 
of the materials in a supercell. This has direct consequences on the nature of the 
magnetic states of superlattices with different thicknesses. 
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1 Introduction 



The magnetic properties of magnetic thin films and multilayers have attracted 
increased attention in recent years [1]. This is because of the easiness of their 
preparation made possible by the recent advances of modern vacuum science, in 
particular the molecular beam epitaxy technique, the novel magnetic phenom- 
ena they exhibit compared to bulk materials, and their potential technological 
importance. The unusual magnetic properties of these artificial materials are 
mostly attributed to surface and interface effects, to the reduced dimensionality 
and to magnetic interactions. 

The thickness dependence of most of the properties of these materials is now 
an established fact. For instance, magnetization measurements of a series of 
Fe films [2] show that the magnetic moment of Fe depends on the thickness of 
the film. In ultrathin Ni films grown on Cu(OOl) [3], a sharp spin reorientation 
transition at small film thicknesses (from 5 to 7 monolayers) from in-plane to 
perpendicular magnetization was observed. At much higher thicknesses (from 35 
to 70 monolayers) the magnetic moments reorient gradually and become again 
parallel to the surface. In superlattices with magnetic layers separated by non 
magnetic spacer layers, the coupling between two magnetic layers oscillate gen- 
erally in sign as a function of the thickness of the spacer. This is the case, for 
example, of the epitaxially grown Gd/Y superlattices [4]. Also, a systematic 
study of EuTe/PbTe superlattices [5], shows considerable changes of the Nel 
temperature depending on the number of monolayers of EuTe and PbTe in a 
supercell. 

The experimental activity on these new materials has been accompanied by 
significant theoretical effort in order to understand their novel properties. Be- 
sides, theoretical work has been devoted to the investigation, on thin film and 
multiplayer geometries, of known interesting models. This has led to the predic- 
tion of new interesting physics. 

The Blume-Capel (BC) model for spin 1 is an Ising system which takes into 
account the single-ion crystal field anisotropy. It has been investigated in detail 
using many approximate methods, namely mean field approximation [6,7], high 
temperature series expansion [8], constant-coupling approximation [9], Monte- 
Carlo [10] and renormalization-group [11] techniques. All of these approximate 
schemes suggest the existence of a tricritical point in the phase diagram of the 
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model and the instability of the ordered phase beyond a critical value of the 
anisotropy parameter. On the other hand the BC model has been investigated 
on semi-infinite lattices with modified surface coupling. Several approximate 
methods have been used [12] (and references therein). All of them show the 
possibility to have a phase with ordered surface and disordered bulk, and show 
the existence of the so-called 'extraordinary', 'surface' and 'ordinary' transitions, 
and of multicritical points called special points [13]. More recently, with the 
increased interest in thin films and multilayers, many works have been devoted 
to the influence of the crystal field on the phase transitions of magnetic films 
[14-15] and superlattices [ [16-17]. 

In an earlier paper [18], we considered an Ising spin-1 multilayer (having L 
layers) with an alternating crystal field: Ai on the odd number layers and A2 
on the even number ones. A rich phase transition behavior has been reported 
there and some novel properties( such as the absence of the tricritical point in the 
paramagnetic-ferromagnetic transition line and the stability of the ferromagnetic 
phase at low temperatures for some regions of the parameters of the model, the 
appearance of the reentrant and double reentrant phenomena and the presence 
of multicritical and critical end-points) have been found. In the present work, we 
address our attention to the physical question of the thickness dependence of all 
those results. The system we consider is a superlattice consisting of two different 
ferromagnetic materials A and B of different crystal fields: Ai in A and A2 in 
B. The supercell contains two slabs of materials A and B having different thick- 
nesses: M and N layers respectively. To get rid of complications due to effects 
which we don't want to take into consideration (boundary conditions at the sur- 
face layers and the parity of the total number of the slabs of the two materials) 
we consider a three-dimensional superlattice. We find that many of the results 
mentioned above change qualitatively and /or quantitatively with the variation 
of the thicknesses of the slabs in a supercell, and that novel phenomena may take 
place. The rest of the paper is organized as follows. In section 2 we present the 
model and give the mean field equations of the different order parameters. In 
section 3 we present the results. Section 4 is devoted to a conclusion. 
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2 Model and method 



The system we consider (Fig.l) is a superlattice consisting of two different fer- 
romagnetic materials A and B stacked alternately: the unit cell consists of an 
arbitrary number M of layers of matarial A with the crystal field Ai, and an 
arbitrary number N of layers of material B with crystal field A2 . For simplicity, 
we focus our attention on the simple cubic structure. The Hamiltonian of this 
system is given by 

where Si = ±1, 0. The first sum runs over all pairs of nearest-neighbors and Aj is 
the crystal field which takes the value Ai on the layers of A and A2 on the layers 
of B. The periodic condition suggests that we have to consider only the unit cell 
of thickness L=M+N. The mean field equations of state are straightforwardly 
obtained and are given by 

2sinh/3% 

mk = {Sk) = e^ + 2cosh^ fc (2) 

2 2cosh/3/i fe 
qk = {Sk) = e /3A fc+2 cosh/3/ lfc (3) 

with 

h k = J (m k -i + zm k + m k+ i), for k = 1, M + N (4) 

and 

hk+M+N = h k (5) 

to ensure the periodic condition 

m k+M +N = m k (6) 

In these equations, m k (q k ) is the reduced magnetization (quadrupolar moment) 
of the k th layer, z is the interlayer coordination number (z=4 for the case of 
the square layers considered here). The reduced total magnetization is m = 
jf J2k=i m k- The reduced free energy of the film is given by 

^ M+N ^ M+N 

f = -- ln 1 + 2e _/3Afc cosh/3/i fc + - ^ ( m k-i + zm k + m k+1 )m k 



(7) 
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3 Results and discussions 



We will be interested in the phase diagrams in the (t, d 2 ) plane for different 
ranges of variation of d\ (d\ = Ai/J, d 2 = A 2 /J and t = ksT/J are the 
values of the reduced crystal fields and temperature respectively). In order to 
determine exhaustively the different phase diagrams, we first determine exactly 
the ground state phase diagram in the (d±,d 2 ) plane by looking for the lowest 
energy configurations. For arbitrary numbers of layers M and N of the slabs of 
the two materials A and B in the unit cell, we distinguish four regions of this 
plane (Fig.2): 

(a) di < 3 - 1/M, d 2 < 3 - l/N and M(di - 3) + N(d 2 - 3) < 0. The stable 
configuration is 1 M 1 N , where the spins of the sites of the M layers of the 
slab A and the N layers of the slab B are all equal to 1. 

(b) di < 3 - 1/M, d 2 > 3 - l/N. The stable configuration is l M N , where the 
spins of the slab A are all equal to 1 while those of the slab B are all equal 
to 0. 

(c) di > 3 - 1/M and d 2 < 3 - l/N. The stable configuration is M 1 N . 

(d) di > 3 - 1/M, d 2 > 3 - l/N and M(d 1 - 3) + N(d 2 - 3) > 0. The stable 
configuration is 0*^0^. 

One of the specificities of the superlattice geometry is worth to note here. 
Material A in bulk with d\ > 3 is non magnetic at zero temperature, but inside 
a superlattice it may be ferromagnetic up to d\ = 3 + 1/M (Fig.2). This is an 
interaction and reduced dimensionality effect: d\ may take higher and higher 
values for thinner and thinner slabs of material A. 

From Fig.2 it is clear that for fixed values of d\, the system exhibits various 
transitions by varying d 2 . For the following ranges of variation of d\ : i) < d\ < 
3 - 1/M, ii) 3 - 1/M < di < 3 + 1/M, iii)cq > 3 + 1/M, the system presents the 
following phase transitions respectively: 

i) from the state 1^1^ to the state 1 M N at d 2 = 3 + l/N. 

ii) from the state 1^1^ to the state M N at N(d 2 - 3) = -M (di - 3). 

iii) from the state 0*^1^ to the state M N at d 2 = 3 - l/N. 
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The topology of the ground state phase diagram is the same for arbitrary 
values of M and N (Fig.2). Nevertheless, the regions defining the different ground 
states in the (di,^) plane are M and N dependant. This implies possible phase 
transitions, at zero temperature, due to thickness variations. For example, for 
d\ = 2.6 and ofo = 2.85, the system with (M = 1, N = 2) is non magnetic. With 
increasing the thickness of the layers of material A, the system ((M = 2, N = 2) 
for example) becomes ferromagnetic. This is formally due to the variation of the 
slope of the line (oblic segment of equation M(d\ — 3)+N(d2 — 3) = 0) separating 
the ferromagnetic from the the non magnetic states, when M or N change. 

For finite temperatures the transition lines are determined by solving numer- 
ically Eq.(2) which may have more than one solution. The one which minimizes 
the free energy (Eq.(7)) corresponds to the stable phase. Transitions of the 
second order are characterized by a continuous vanishing of the magnetizations 
while those of the first order exhibit discontinuities at the transition points. 

In order to describe the different entities in the phase diagrams, the Griffiths 
notations [19] will be adopted: the critical end-point B m A n denotes the intersec- 
tion of m lines of second order and n of first order. The multicritical point B m 
denotes the intersection of m lines of second order. The tricritical point, which 
is the intersection of a line of second order and a line of first order, is denoted in 
particular by C. 

In addition, each phase will be represented by the juxtaposition of two ex- 
pressions each one giving the state of a slab A or B in the unit cell. The notation 
(F M )(FP0 N ~ 2 PFP) means that the slab A is in the state (F M ) (i.e. its M lay- 
ers are all in ferromagnetic states), and the slab B is in the state (F p N ~ 2p F p ) 
(i.e. the first p layers of B, and the last p ones because of symmetry, are all 
in ferromagnetic states while the N — 2p inner ones are all in the non magnetic 
state 0) . A layer in the non- magnetic state has a vanishing magnetization and 
quadrupolar moment. In the paramagnetic state P only the quadrupolar mo- 
ment is non vanishing, and in the ferromagnetic state F both quantities are non 
vanishing. 

In what follows, we will give the phase diagrams in the (t,^) plane for given 
values of M and N (M = 5, N = 9). Meanwhile we give the phase diagrams for 
other values of M or N if their topologies are different from those of (M = 5, 
N = 9). For arbitrary fixed values of M and N, the different phase diagrams 
may be gathered into four groups, each one is characterized by a certain range 
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of variation of d\ : 



1) For < di < d trc (e.g. Fig.3(a)), the system presents the paramag- 
netic phase (P M )(P N ) and the ferromagnetic ones (F M )(F N ) (called the or- 
dered phase) and (F M )(F p N ~ 2p F p ) (called the partially ordered phases). The 
system transits from one ferromagnetic phase to another when a certain layer, 
and its symmetric, transits from the non magnetic state to a ferromagnetic one 
F. There is a critical line separating the paramagnetic phase (P M )(P N ) from 
the ferromagnetic ones. The ordered phase is separated from the partially or- 
dered ones by a transition line of first order at low temperatures, which meets 
a critical line at a tricritical point. The critical line meets the paramagnetic- 
ferromagnetic (P-F) transition line at a multicritical point of the type B s . The 
partially ordered phases are separated by critical lines, each one representing 
the transition of a layer, and its symmetric, from a non magnetic to a magnetic 
state. These lines end, at low temperatures, at some tricritical points and end 
at some multicritical points of the type -B 3 on the P-F transition line. They 
furthermore present the reentrant behavior near the P-F transition. Note that 
the first order line, which meets the d^ axis at d2=3+l/N, presents an inclina- 
tion inside the ordered phase (F M )(F N ). This means that, contrary to the zero 
temperature case, some layers of the slab B may be non magnetic even for cfo 
less, but close to 3+1/N, at low temperatures. This also means the existence of 
the reentrant behavior where for d2 close to 3+1/N the system, at the beginning 
in the ordered phase (F M )(F N ), transits to some partially ordered phases and 
then returns to the ordered one. This is clarified by the inset in Fig. 3(a). It has 
been found that this behavior does not exist for small values of N (e.g. N=2) 
in which case the first order line mentioned above is parallel to the t-axis. The 
topology described here is characterized, contrary to the ordinary Blume-Capel 
model, by the absence of a tricritical point in the (P-F) transition line. The value 
dfr^ of d\ above which this point appears is found to be the mean field value 
of the tricritical crystal field of the spin-1 Blume-Capel film with M layers and 
a uniform crystal field. d[^} increases from d[l\ = 2. for one layer to d^) = 3. 
for the bulk. The M dependance of d[^} implies that for particular values of 
d\, superlattices with M and M' so that d^} < d\ < d[^ c \ have qualitatively 
different phase diagrams. For example, for d± = 2.62, the system with thickness 
M = 4 (dlfl = 2.59) presents a tricritical point, but not the system with M = 5 
(dfr) c = 2.65). The absence of the tricritical point in the phase diagrams of the 
Ising mixed spin — 1/2, spin — 1 superlattice of [17], may be considered as the 
particular case d\ = 0. of the present model. 



7 



2) For d\Z <d 1 < 3-1/ M, three types of topology are found with increasing 
di, all of them present a tricritical point in the P-F transition line and present the 
same features of the type-1 topology with some differences. For d\ in the vicinity 
of dpfj some of the lines of transition (F-F lines) separating the ferromagnetic 
phases end at some multicritical points of the type B 3 on the (P-F) transition 
line, but some others end at points of the type BA 2 [e.g. Fig.3(b)]. For higher 
values of d\, all the F-F lines end at the multicritical lines of the type BA 2 on the 
P-F line [e.g. Fig. 3(c)]. Finally for d\ close to (3-1/M), some partially ordered 
phases disappear ((i ?5 )(i ?4 0i ?4 ) in the case of M=5, N=9). This is true only for 
thick slabs of material B and means that the inner layers of this material can not 
be driven to the ferromagnetic state by the temperature. In the other hand, the 
transition line separating the ordered and the partially ordered phases is now of 
totally first order, and meets the (P-F) line at a triple-point A 3 [e.g. Fig. 3(d)]. 
Such topology does not exist for small values of N (N < 7). 

For d± < 3 — 1/M the ferromagnetic phase is stable even for an infinitely 
large value of di at temperatures lower than a temperature t c ^ M \d\) which is 
independent of the thickness N of the slab B. t c ^ M \d\) is but the mean-field 
critical temperature of the spin-1 Blume-Capel film with M layers and a uniform 
crystal field that is equal to d\. This may be understood beyond the mean-field 
approximation [13]. tJ- M \d\) is a decreasing function of d\. This means that for 
large values of di the (P-F) transition line takes place at lower temperatures for 
higher values of d\. In fact, as is clear from Figs.2(a-d), this is true for all di- 
This fact explains the progressively disappearance of the re-entrance in the lines 
separating the ferromagnetic phases near the (P-F) transition line, the eventual 
disappearance of some partially ordered phases, and the eventual appearance of 
the A 3 point. 

For d\ > 3 — 1/M, the ferromagnetic phase is unstable for large values of 
di-, as will be seen below, in accordance with the ground-state phase diagram. 
This implies that the magnetic state at low temperatures is M dependant. For 
example, for d\ = 2.75, the system with (M = 5, N = 4) exhibits a ferromagnetic 
order at low temperatures for arbitrary values of di, but not the system with 
(M = 2,N = 4). 

3) For 3 — 1/M < d\ < 3 + 1/M, no layers transitions may take place. The 
phase diagram for (M = 5, N = 9) is of the same type as that of the ordinary 
spin-1 BC model [e.g. Fig. 3(e)]. The first order line reaches the di axis at 
di given by N(di — 3) = — M(d\ — 3). Nevertheless, a reentrants and double 
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reentrants behaviors may appear in the P-F transition line for varied values of 
the thicknesses. This is the case for values of d\ close to (3 — 1/M) for small 
values of M relatively to N (e.g. Fig.4(a) for M=2 and N=9), and for values of 
di close to (3 + 1/M) for high values of M relatively to N (e.g. Fig.4(b) for M=9 
and N=2). 

4) For d\ > 3 + 1/M many types of topology are found for M=5 and N=9, all 
of them have a (P-F) transition line containing a tricritical point C and meet the 
^2 axis at o?2 = 3 — 1/N. Some of the ferromagnetic phases (0 5 )(F 9 ), (F0 3 F)(F 9 ), 
(F 2 0F 2 )(F 9 ) or (0 5 )(-F 9 ) may appear in the phase diagrams depending on the 
values of d±. For d± close to (3+ 1/M), all of these phases exist and are separated 
by critical lines of approximately constant temperatures which meet the P-F 
transition line at three critical end points of the type (BA 2 ) [e.g. Fig. 3(f)]. With 
increasing di, the following changes take place, 

-the ferromagnetic phase (i ?5 )(i ?9 ) becomes totally inside the phase (F 2 0F 2 )(F 9 ) 
which exhibits an island behavior [e.g. Fig.3(g)]. 

- The phase (F 2 0F 2 )(F 9 ) becomes itself totally inside the phase (F0 3 F)(F 9 ) 
and two islands then appear [e.g. Fig.3(h)]. 

-The ferromagnetic phase (F 5 )(F 9 ) disappears from the phase diagram [e.g. 
Fig.3(i)] 

- The phase (F 2 0F 2 )(F 9 ) disappears too, and only the two phases (0 5 )(F 9 ) 
and (F0 3 F)(F 9 ) still exist. They are separated by a critical line which meets 
the (P-F) transition line at a (BA 2 ) point [e.g. Fig.3(j)]. 

- The latter line exhibits a reentrant behavior and the phase (F0 3 F)(F 9 ) 
becomes totally inside the (0 5 )(F 9 ) phase [e.g. Fig.3(k)]. 

- The phase (F0 3 F)(F 9 ) disappears and the topology becomes of the ordinary 
BC type, but now with the partially ordered phase (0 5 )(-F 9 ) at low temperatures, 
instead of the ordered phase (F 5 )(F 9 ) [e.g. Fig.3(l)]. 

For general values of M and N, analogous phase diagrams are found, the 
number of the ferromagnetic phases which appear in these diagrams depends 
obviously on M and N and depends on d\. For some values of M and N, the 
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(P-F) transition line may exhibit the reentrant or double-reentrant behaviors for 
di close to 3+1/M. This is the case, for example, for M=9 and N=2 (Fig.4(b)). 



4 Conclusion 



We have studied, using the mean field theory, the phase diagrams of a spin- 
1 Blume-Capel superlattice consisting of two different ferromagnetic materials 
with two different crystal fields Ai and A2. The ground state phase diagram in 
the (di, d%) plane was determined exactly for arbitrary values of the thicknesses M 
and N of the constituents in a unit cell. At finite temperature, the phase diagrams 
have been given in the (t,^) plane for different ranges of variation of d\ .They 
exhibit a variety of multicritical points, in particular the usual tricritical point C 
in the paramagnetic-ferromagnetic line of transition. This point does not appear 
for values of d\ below a threshold value d[^} , which depends only on the thickness 
M of the constituent with crystal field Ai in the unit cell. Moreover, lines of 
transition presenting the reentrant, double reentrant and island behaviors may 
also take place. The phase diagrams depend qualitatively and/or quantitatively 
on M and N. This has direct consequences on the magnetic nature of superlattices 
with different thicknesses. 
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Figure captions 

Fig.l A schematic illustration of the model. 

Fig.2 The ground state phase diagram in the (d^cfe) plane for arbitrary num- 
bers M and N of the layers of the slabs A and B respectively, in a unit 
cell. 

Fig. 3 The phase diagrams in the (t,^) plane for M = 5 and A = 9 and generic 
values of di. (a) di = 1., (b) di = 2.66, (c) di = 2.7, (d) d-y = 2.79, (e) 
di = 3., (f) d 1 = 3.5, (g) dx = 5., (h) d x = 6., (i) di = 7., (j) di = 8.7, (k) 
di = 13., (1) di = 18. The solid lines are of second-order, the doted ones 
are of first-order. Filled circles are the tricritical points, empty circles are 
the -B 3 multicritical points, filled squarres are the BA 2 critical end-points 
and the empty squarre is the A 3 end-point. 

Fig.4 The phase diagrams in the (t,d2) plane for (a) M = 2 and N = 9 
(3 - 1/M = 2.5), (b)M = 9 and N = 2 (3 + l/M = 3.11) .The number 
accompanying each curve is the value of d±. The solid lines are of second- 
order, the doted ones are of first-order. Filled circles are the tricritical 
points. 
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